Entanglement Spectrum and Entanglement Thermodynamics of Quantum Hall 

Bilayers at u = 1 
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We study the entanglement spectra of bilayer quantum Hall systems at total filling factor f — 1. 
In the interlayer-coherent phase at layer separations smaller than a critical value, the entanglement 
spectra show a striking similarity to the energy spectra of the corresponding monolayer systems 
around half filling. The transition to the incoherent phase can be followed in terms of low-lying 
entanglement levels. Finally, we describe the connection between those two types of spectra in 
terms of an effective temperature leading to relations for the entanglement entropy which are in full 
analogy to canonical thermodynamics. 

PACS numbers: 73.43.-f,73.43.Nq 



I. INTRODUCTION 

Experimental studies of quantum Hall bilayers at to- 
tal filling factor 1^ = 1 have revealed a series of in- 
triguing phenomena-"^. From a theoretical point of 
view, such strongly interacting systems pose intricate 
and subtle problems of many-body physics, and many 
unexpected observations are still lacking a fundamen- 
tal understanding!. Here the exact numerical treatment 
of small systems has provided important insights and 
guidance of theoretical intuition^rJ^. Moreover, in the 
last decade many-body physics has been enriched by the 
concept of entanglement, whose extensive and system- 
atic study originated in the field of quantum information 
theory^^. In particular, the notion of the entanglement 
spectrum has led to novel insights in the physics of quan- 
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turn Hall monolayers at fractional filling factor; 
systems of one^ir— and twc^*' spatial dimensions, quan- 
tum hall insulators2!i2^, and also rotating Bose-Einstein 
condensates^^. The entanglement spectrum of a bipartite 
system is defined in terms of the Schmidt decomposition 
of its ground state |?/;), 



(1) 



where the states ( form an orthonormal basis 

of the subsystem A (B), and the nonnegative quantities 
are the dimensionless levels of the entanglement spec- 
trum. In the present paper we study the entanglement 
spectrum of quantum Hall bilayers at total filling factor 
!/ = 1, where the subsystems are naturally defined by 
the two layers. Our investigations are based on exact 
numerical diagonalizations in the spherical geometry^S. 



II. RESULTS AND DISCUSSION 

Depending on the layer separation and the ampli- 
tude for particle tunneling between the layers, quantum 
Hall bilayers exhibit a quantum phase transition between 
a phase showing spontaneous interlayer coherence and 



quantized Hall conductance, and a phase where such ef- 
fects are absent with the two layers being merely un- 
correlated. The corresponding phase diagram was first 
mapped out experimentally by Murphy et al} and later 
reproduced theoretically via numerical diagonalizations^. 
The latter study relied on a careful analysis of fluctu- 
ations of the pseudospin describing the layer degree of 
freedom. Moreover, in the case of of very small tunnel- 
ing amplitude, the phase transition is also signaled in 
experiment by the generation of a pronounced peak in 
the tunneling conductance at zero bias voltage^. Here 
the critical value of the layer separation of d = 1.83^ 
[t. magnetic length) was nicely reproduced by numeri- 
cal calculations taking into account the finite width of 
the two wells^ii^i^. In the following, however, we shall 
for simplicity concentrate on the case of zero well width 
where the transition is found numerically at d = 1.3^, 
and we will focus on the situation of vanishing tunnel- 
ing amplitude. Moreover, we will assume the physical 
electron spin (as opposed to its pseudospin) to be fully 
polarized by the perpendicular magnetic field such that 
it does not play a role in the present study. 

In the limit of zero layer separation d = both layers 
merge to a fully filled monolayer whose Coulomb ground 
state is given by the well-known Slater determinant of 
the Vandermonde type, 
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where the Zi are the complex electron coordinates in 
the planar geometry. An analogous wave function ex- 
ists in the spherical model of N electrons where the 
sphere is penetrated by — 1 flux quanta^S. Here 
the z-component of the pseudospin can take values of 
= —N/2, . . .N/2 meaning that, say, the top/bottom 
layer contains N/2 ± electrons. The case = cor- 
responds to Halperin's 111 wave function'^-. Tracing out 
one of the layers means summing over N /2 ± electron 
coordinates in the pure state p = \'ip){ip\ with being 
a single Slater determinant such that the resulting N/2- 
body reduced density matrix is proportional to the unit 
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TABLE I: Angular momentum !/„ of the low-lying entangle- 
ment levels at various system sizes A'^ of the underlying 
bilayer system. 



matrix with all eigenvalues being equal to 1/(7v/2±t^)- 
Thus, at d = the entanglement spectrum shrinks to a 
single value. 

Let us now turn to the case of finite c? > 0. The 
Coulomb ground state of of the bilayer system in spher- 
ical geometry has vanishing total angular momentum 
i = 0. Therefore, the reduced density matrix is also 
rotationally invariant, and its eigenvalues occur in multi- 
plets of the angular momentum. Fig. [1] presents entangle- 
ment spectra obtained from balanced bilayers (T^ = 0) 
at small but finite layer separation d = 0.1£ for vari- 
ous system sizes A'^ and compares them with the energy 
spectra of half-filled monolayers containing Nm = N/2 
particles. As seen from the figure, both types of spectra 
show a striking similarity which is, on the total scale of 
the spectra, particularly evident at smaller system sizes. 
This observation also continues to larger layer separa- 
tions as shown in Fig. [5] for the cases d — 0.71 and 
d = With increasing layer separation the entangle- 
ment spectra widen up (while shrinking to a single level 
for d — ^ 0), but remain similar in shape to the mono- 
layer energy spectra. At layer separations exceeding the 
critical value d = l.M a qualitatively different behav- 
ior sets in: The global shape of the entanglement spec- 
trum remains reasonably close to the monolayer energy 
spectrum, but differences are developed in details of the 
data. This can be seen most clearly by following the low- 
est entanglement level and the three following levels 
£.71, n e {1, 2, 3}, as a function of the layer separation'^2,. 
Those entanglement levels (along with their counterparts 
in the monolayer energy spectra) occur at different values 
of the angular momentum L„ which depend on system 
size and are summarized in table HI Fig. [3] shows these 
low-lying entanglement levels as a function of the layer 
separation. While at d < all levels decrease with in- 
creasing layer separation, the properties of the spectrum 
change qualitatively at d « 1.3^: Here the excited levels 
start to increase with the layer separation while the low- 
est level approaches ~ ln(2Lo -I- 1). Indeed, at large 
layer separations the reduced density matrix of a single 
layer is dominated by a multiplet with total angular mo- 
mentum Lq with all other eigenvalues being exponentially 
close to zero. In this sense, the ground state of the bilayer 
system carrying angular momentum quantum numbers 
L = = can at large layer separation be viewed to be 
composed of two monolayers with L = Lq according to a 
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leading to a reduced density matrix of the form 
\ i 

Pred^ y] — --r\Lo,m){Lo,m\. 
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Thus, the entanglement entropy S = {—hipped), (■) '■= 
tr(-/9red), interpolates, similarly as the lowest entangle- 
ment level ^0, between a value of 5 = In (^2) at ^ = 
and S — ln(2Lo + 1) at large layer separation. In the left 
bottom panel of Fig. 0] we have plotted S along with its 
variance 
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as a function of d for a bilayer system of = 12 electrons. 
The entropy S shows an inflection point, accompanied by 
a maximum of AS*, near d — lAi, which is close to the 
finite-size value of the phase boundary obtained earlier 
from an analysis of the pseudospin fluctuations^!^. 

The observation that the entanglement spectrum of 
interlayer-coherent quantum Hall bilayers shows an in- 
triguing similarity to the the energy spectrum of the cor- 
responding monolayer extends also to unbalanced sys- 
tems. This is illustrated in the top panels of Fig. |4] 
where the entanglement spectra of A^ = 12 electrons 
with e {1,2} obtained by tracing out the top layer 
are compared with the corresponding monolayer spec- 
tra. Entanglement spectra obtained from bilayer sys- 
tems with negative (or, alternatively, by tracing out 
the bottom layer) are related to the previous ones by a 
particle-hole transformation are therefore identical, and 
the corresponding monolayer energy spectra just differ 
by an additive constant. 

The investigations so far have concentrated on the case 
of vanishing single-particle tunneling between the lay- 
ers. At finite tunneling amplitude, the particle number of 
each layer is no longer a good quantum number, and the 
monolayer Fock space consists of invariant subspaces of 
the reduced density matrix which are characterized both 
by total angular momentum and particle number. There- 
fore, for not too large tunneling, the entanglement spec- 
tra within each subspace of given particle number will 
be similar to those given in Figs I1|2I3I up to an additive 
constant describing the spectral weight of the respective 
subspace. 

Focusing again on the case of zero tunneling, our re- 
sults suggest that the reduced density matrix fulfills, for 
not too large layer separation d, the following approxi- 
mate relation, 
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where Z = tr exp(— /^Hm), and Hm describes the 
Coulomb repulsion in a half-filled monolayer with A^ flux 
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quanta and N„i = N/2 electrons. The (inverse) "entan- 
glement temperature" /3 is a parameter depending on d, 
and the occurrence of the partition function Z in Eq. ^ 
ensures the condition trp^ed = 1- Clearly, d ^ im- 
plies (3^0 and Z = (^2)- ^^'^ bottom panel of 
Fig.|3]shows least-mean-square results for /3 as a function 
of layer separation. As seen there, the data depend only 
very weakly on system size. Moreover, defining a "en- 
tanglement free energy" by F = — TlnZ, T = 1/(3, one 
derives from Eq. ^ the familiar thermodynamic relation 
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with E = {Hm)- Thus, the entanglement entropy of the 
bilayer system turns out to have an immediate thermo- 
dynamic meaning for the monolayer system characterized 
by the parameter T = These findings suggest to de- 
fine, in full analogy with canonical thermodynamics, a 
"specific heat" via 



r^_rj.dS_ dS 
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In the mid bottom panel of Fig. U we have plotted C as 
a function of layer separation for a system of = 12 
electrons, as seen, the specific heat has a pronounced 
peak near the phase boundary which should be seen as a 
strong hint for the underlying phase transition being first 
order. Finally, from Eq. ^ one also derives the result 



(9) 



which provides a direct link between the variance of 
the entanglement entropy of the bilayer system and the 
energy variance of the monolayer. As seen in Fig. Ul 
AiS* shows a maximum at the phase boundary which 
should be considered as another piece of evidence that, in 
the infinite- volume limit, the phase transition is of first 
order— rather than continuousii, i.e. the entropy changes 
discontinuously. 

Very recently, Poilblanc has performed a study of en- 
tanglement spectra of Heisenberg spin ladders, where a 
similarly striking resemblance to the energy spectra of 
single spin chains was found^^. This observation could 
also be described in terms of an effective temperature. 
The results of both papers raise the question of a com- 
mon underlying mechanism being at work, also extend- 
ing previous studies on fractionally filled quantum Hall 
monolayersJi^^° to bilayers at filling factor v = 1 inves- 
tigated here. Most interestingly, the latter direction of 
work makes immediate contact to the enigmatic system 
of a half-filled quantum Hall monolayer. However, here 
we should add the caveat that the monolayer systems 



that naturally arise in our present work contain N elec- 
trons in = 2N — 1 flux quanta whereas in many previ- 
ous numerical studies on (electron-spin polarized) mono- 
layers at half filling considered using spherical geometry 
other relations of the form A^^ = 2N — M, M being a 
small integer, have been considered^^. To eliminate such 
a finite-size "shift problem" it would be desirable to re- 
peat the calculations reported on here in the toroidal 
geometry. 

We note that a different recipe of how to eliminate 
finite-size properties from entanglement spectra of quan- 
tum Hall monolayers in the spherical geometry was put 
forward very recently by Thomale et alJ^. These authors 
propose to strip all normalization factors from single- 
particle states entering the Slater determinants of many- 
body wave functions. This procedure of achieving a 
"conformal limit" conserves azimuthal symmetry on the 
sphere (as appropriate for the situation of Ref J^) , but de- 
stroys full rotational covariance, which is desired in the 
present study on bilayer systems. Thus, the prescription 
proposed in RefJ^ is unfortunately not applicable to our 
investigations here. 



III. CONCLUSIONS 

In summary, we have studied the entanglement spec- 
tra of bilayer quantum Hall systems at total filling factor 
v — 1. Our investigations are based on exact numeri- 
cal diagonalizations using the spherical geometry. In the 
interlayer-coherent phase at small layer separations, the 
entanglement spectra show a striking similarity to the 
energy spectra of the corresponding monolayer systems 
around half filling. The transition to the incoherent phase 
can be followed in terms of low-lying entanglement levels, 
constituting a link between the entanglement spectrum 
and a quantum phase transition. Clear signatures of the 
quantum phase transition are also shown by the entangle- 
ment entropy along with its fluctuation. Moreover, the 
connection between those two types of spectra can be de- 
scribed in terms of an effective temperature which gives 
rise to relations for the entanglement entropy being fully 
analogous to canonical thermodynamics. In particular, 
the specific heat derived from this formalism provides a 
strong hint for the phase transition being of first order—. 
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FIG. 1: Left: Entanglement spectra obtained from balanced 
bilayers at layer separation d = 0.1.^! for various system sizes 
A^. Right: Energy spectra of half-filled monolayers with A'^/2 
electrons. 




FIG. 2: Entanglement spectra of balanced bilayers at layer 
separation d = 0.7^ (left) and d = 1.3^ (right). 
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FIG. 3: The four lowest entanglement levels as a function of 

the layer separation for different system sizes. At the phase 
boundary at d « 1.3£ the entanglement spectrum qualita- 
tively changes. 




FIG. 4: Top panels; Entanglement spectra of = 12 elec- 
trons with G {li2} are compared with the corresponding 
monolayer spectra . Left bottom panel: Entanglement en- 
tropy S along with its variance AS as a function of the layer 
separation d for a balanced bilayer system of iV = 12 elec- 
trons. Right bottom panels: The "inverse temperature" /3 
and the "specific heat" C7 as a function of layer separation. 



